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The fluctuation properties of the human electroencephalo- 
gram (EEG) time series are studied using detrended fluctu- 
ation analysis. For all 128 channels in each of 18 subjects 
studied, it is found that the standard deviation of the fluctu- 
ations exhibits scaling behaviors in two regions. Topographi- 
cal plots of the scaling exponents reveal the spatial structure 
of the nonlinear electrical activities recorded on the scalp. 
Moment analyses are performed to extract the gross features 
of all the scaling exponents. The correlation between the two 
scaling exponents in each channel is also examined. It is found 
that two indices can characterize the overall properties of the 
fluctuation behaviors of the brain dynamics for every subject 
and that they vary widely across the subjects. 

FAGS numbers: 05.45. Tp, 87.19.La, 87.90. -)-y 



I. INTRODUCTION 

The scalp electroencephalogram (EEG) provides a 
wealth of information about human brain dynamics. The 
complex nature of brain dynamics results in a high degree 
of fluctuations in both the spatial and temporal aspects 
of the EEG signals. To extract the salient properties 
from the data is the primary objective of any method of 
analysis. We present in this paper a novel method that 
explores the scaling behavior of the fluctuations and uses 
moment analysis to reduce the complexity of the results 
obtained. 

The most common methods of EEG time series analy- 
ses are event-related time ensemble averaging and Fourier 
decomposition, both of which are based implicitly on 
assumptions of linearity Since the physiological 

mechanisms underlying the scalp EEG are generally non- 
linear, they can generate fluctuations that are not best 
described by linear decomposition. Moreover, the rest- 
ing EEG always displays a broad-banded power spec- 
trum, so in Fourier analysis one must arbitrarily de- 
flne frequency bands (6, 9,a, ■ ■ ■) which may not actu- 
ally delineate different dynamical mechanisms. Wavelet 
analyses have also been applied to examine EEG time 
series but at a sacrifice of the ability to describe 
long-range temporal correlations. Chaos analyses have 
been applied to quantify the nonlinear behavior of the 



brain [Q-^, but typically require a long period of time 
to compute attractor properties for a single time series. 
Moreover, chaos-based approaches assume the existence 
of low-dimensional attractors, and this is probably not 
generally a valid assumption for the brain. In this paper, 
we discuss a method that analyzes the fluctuations in 
EEG over a short period of time (around 10s), and makes 
use of the information conveyed by all 128 channels on 
the scalp. We show the existence of scaling behaviors of 
certain measures of fluctuations in all channels and in all 
subjects. We propose two global measures of the spatio- 
temporal signals that characterize the distinctive nature 
of EEG for each subject. 

The study of scaling behavior emphasizes the relation- 
ship across time scales. We aim to find what is universal 
among all channels as well as what varies among them. 
The former is obviously important by virtue of its uni- 
versality for a given subject; how that universal quantity 
varies from subject to subject is clearly interesting. What 
varies from channel to channel is perhaps even more in- 
teresting, since it has implications for describing focal 
features which could have functional or clinical relevance. 

Our procedure is to focus initially on one channel at a 
time. Thus it is a study of the local temporal behavior 
and the determination of a few parameters (scaling expo- 
nents) that effectively summarize the fluctuation proper- 
ties of the time series. The second phase of our procedure 
is to describe the global behavior of all channels and to 
arrive at two numbers that summarize the variability of 
these temporal measures across the entire scalp surface. 
This dramatic data reduction necessarily trades detail for 
succinctness, but such reduction is exactly what is needed 
to allow easy discrimination between brain states. 

The emphasis in this paper is on the method of analy- 
sis. It is not our aim here to perform clinical and cogni- 
tive analyses. Due to the fact that the EEG data avail- 
able to us are short in time duration and few in the 
number of subjects studied, it is not feasible for us to 
make reliable inference on the physiological implications 
of our flndings. Nevertheless, the data are sufficient for 
the extraction of interesting behaviors that are channel 
dependent as well as subject dependent. 
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II. DETRENDED FLUCTUATION ANALYSIS 

The specific method we use in the first phase is de- 
trended fluctuation analysis (DFA). This analysis is not 
new. It was proposed for the investigation of correlation 
properties in non-stationary time series and applied to 
the studies of heartbeat and DNA nucleotides |^ . It 
has also been applied to EEG but with somewhat 
different emphases than those presented here. Since the 
analysis considers only the fluctuations from the local lin- 
ear trends, it is insensitive to spurious correlations intro- 
duced by non-stationary external trends. By examining 
the scaling behavior one can learn about the nature of 
short-range and long-range correlations. 

Let an EEG time series be denoted by y{t), where t 
is discrete time ranging from 1 to T. Divide the entire 
range of t to be investigated into B equal bins, discard- 
ing any remainder, so that each bin has k = floor(r/i?) 
time points. Within each bin, labeled b [h — 1, • • • , S), 
perform a least-square fit of y{t) by a straight line, y},[t), 
i.e., yh{t) = Linear- fit for (6- l)fc < t < bk. That is 
the semi- local trend for the &th bin. Define F^{k) to be 
the variance of the fluctuation y{t) from yi,{t) in the bth 
bin, i.e., 

bk 

= E iyit)-ydt)r (1) 

t=(6-l)fc+l 

It is a measure of the semi-locally detrended fluctuation 
in bin b. The average of F^{k) over all bins is 

F'{k)^^J2Pbik). (2) 

b=l 

F{k) is then the RMS fluctuation from the semi-local 
trends in B bins each having k time points. 

The study of the dependence of F{k) on the bin size k 
is the essence of DFA . If it is a power-law behavior 

F{k) oc k", (3) 

then the scaling exponent a is an indicator of the corre- 
lations of the fluctuations in EEG, which depends on the 
relationship of these fluctuations across time scales. Since 
DFA considers only the fluctuations from the semi-local 
linear trends, it is insensitive to spurious correlations in- 
troduced by non-stationary external trends. This is a 
practical advantage since EEG acquisition systems often 
suffer from slow drifts associated with gradual changes in 
the quality of electrode contact to the skin. The analysis 
also liberates our result from the dependence on the over- 
all magnitude of the voltage y{t) recorded by each probe, 
which is an advantage since overall signal amplitude can 
vary across subjects, presumably due to differences in 
skull conductivity and other factors. 

Resting EEG data were collected for 18 subjects us- 
ing a 128-channel commercial EEG system, with scalp- 
electrode impedences ranging from 10 to 40 k^l. The 



data were hardware filtered between 0.1 and 100 Hz, then 
digitized at 250 points/sec. After acquisition, T « 10s 
lengths of simultaneous time series in all channels are 
chosen, free of artifacts such as eye blink and head move- 
ments. At each time point, the average across all elec- 
trodes was subtracted, to remove approximately the ef- 
fect of the reference electrode We investigate the 
range of k from 3 to 500 in approximately equal steps of 
Infc. 

In Fig. 1 we show three typical time series y{t) in three 
widely separated channels for subject A, labeled 1-3, for 
brevity. While it is clear that both channels 2 and 3 
have substantial 10 Hz oscillations after 0.2s, it is much 
less apparent that there exist any scaling behaviors in 
all three channels. The corresponding values of F{k) are 
shown in the log-log plot in Fig. 2. Evidently, the strik- 
ing feature is that there are two scaling regions with a 
discernible bend when the two slopes in the two regions 
are distinctly different. With rare exceptions this feature 
is found in all channels for all subjects. Admittedly, the 
extents of the scaling regions are not wide, so the behav- 
ior does not meet the qualification for scaling in large 
critical systems or in fractal geometrical objects. How- 
ever, since the behavior is so universal across channels 
and subjects, and since the temporal scales involved are 
physiologically relevant, the scaling behavior is a feature 
of EEG that conveys an important property of the brain 
activity and should not be overlooked. 

III. SCALING AND NONSCALING PROPERTIES 

The fact that there exist two scaling regions suggests 
that the lack of scaling in the region between the two im- 
plies the existence of some significant time scale. From 
Fig. 1 one indeed sees roughly periodic oscillations in Ch. 
2 and 3. One may therefore be tempted to think that if, 
instead of considering the fluctuations from the linear 
semi- local trends yb{t), one studies the fluctuations from 
periodic oscillations, then the bend would disappear and 
the two scaling regions might be joined to become one. 
However, even if that were true, such a procedure should 
not be used for two reasons. First, not all channels ex- 
hibit obvious oscillatory behaviors with definite frequen- 
cies. Channel 1 in Fig. 1 is one such example. Whatever 
detrending one chooses should be universally applied to 
all channels in order to avoid introducing discrepencies 
across the channels due to external intervention. Sec- 
ond, to determine the frequency of the oscillatory trend 
requires a Fourier analysis, which is precisely what our 
approach attempts to circumvent. To decide on a sinu- 
soidal wave of a particular frequency as reference for de- 
trending involves arbitrariness and is unlikely to lead to 
any simplification in the global picture. The simplest and 
unbiased approach is to use the semi-local linear trends, 
as we have done. 

To quantify the scaling behavior, we perform a linear 
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fit in Region I for 1 < Infc < 2.5 and denote the slope by 
ai, and similarly in Region II for 3.5 < Infc < 5.75 with 
slope denoted by a2- Visual inspection for each of the 18 
subjects verifies that fitting this way does a remarkably 
good job of characterizing the slopes in the two regions. 
Knowing the two straight lines in each channel allows us 
to determine the location of their intercept, Iuk, which 
gives a good approximation for the position of the bend 
in Infc. We find that, whereas ai and a2 can fluctuate 
widely from channel to channel, k is limited to a nar- 
row range in most subjects. The average value of ln«; for 
each subject ranges from 2.6 to 3.6, with a grand aver- 
age across subjects to be approximately 3.1. It should, 
however, be noted that when ai and a2 are nearly the 
same, as is the case for Ch. 1 in Fig. 2, the determination 
of K, by the intersection of the two straight lines is not 
reliable. Nevertheless, it is visually clear that the bend 
occurs in the vicinity of Iuk = 3.1. 

Since scaling behavior means that the system exam- 
ined has no intrinsic scale, scale noninvariance at k, im- 
plies that K is related to a dominant frequency of oscilla- 
tion in the time series. It is at this point that a contact 
can be made with the usual Fourier analysis. Although 
our analysis focuses on scale invariant quantities, i.e., the 
dimensionless scaling exponents, it is worth digressing 
momentarily to establish this contact. To do this, we 
loosely associate the time scale k with the period of a 
sine wave with frequency /. If the data acquisition rate 
is denoted by r, then the frequency / corresponding to k 
is 

/ = t/k. (4) 

For our data acquisition we have r = 250 points/sec. For 
the across-subject average of Iuk = 3.1, we get from Eq. 
(4) / = 11.3 Hz. That is in the middle of the traditional 
a (8-13 Hz) EEG frequency band. Thus the dominant 
periodic oscillation apparent in Fig. 1 does reveal itself 
in the study of the scaling behavior. If one's interest is 
in the frequency content of the EEG time series, then 
Fourier analysis is more direct. However, if the interest 
is in the fluctuation properties and their relationships 
across time scales, then DFA is more effective. Hereafter, 
frequency will play no essential role in the remainder of 
this paper. 

For each subject we have 128 pairs of values of (ai, 02), 
which summarize the temporal fluctuations in terms of 
scaling exponents. In Fig. 3, we exhibit by scatter plot 
the values of ai and a2 of all channels for subject A. The 
three points marked with circles correspond to the chan- 
nels shown in Fig. 2. The error bars in Fig. 3 indicate the 
goodness of fit of the two regions by straight lines. Since 
the variability of the scaling exponents across channels 
is large compared to the error in fit, the different val- 
ues convey numerically meaningful information. For this 
subject. Region I exhibits better scaling behavior than 
Region II, although the error bars on a2 are not so large 
as to call into question the power-law description in Re- 
gion II. 



Overall, for subject A, the scaling exponents are in 
the ranges: 0.19 < ai < 1.44 and 0.018 < < 0.489. 
Whereas ai is widely distributed, 02 is sharply peaked 
at 0.1 and has a long tail. The value of a = 0.5 corre- 
sponds to random walk with no correlation among the 
various time points. For a ^ 0.5 there are correlations: 
Region I corresponds to short-range correlation. Region 
II long-range, with n giving a quantitative demarkation 
between the two. In most channels we find ai > 012, al- 
though there are a few where ai ~ a2. The scatter plots 
of all other subjects are similar in general features to the 
one shown in Fig. 3, but vary in detail from subject to 
subject. It is impractical to show them all in this pa- 
per. Evidently, it is desirable to find a way to quantify 
succinctly these 128 pairs of numbers so that one can 
effectively compare the results across subjects. 

A scatter plot such as Fig. 3 reveals very well how the 
ai exponents of all the channels are related to one an- 
other. However, it shows nothing about the locations of 
the channels on the scalp. To show that, we use the to- 
pographical plots of ai and a2 separately, as in Fig. 4, 
to exhibit the spatial structure of the signals extracted. 
The dissociation of ai from 02 is a price paid to gain the 
spatial perspective on the scalp. Topographic plots such 
as this may be useful for specifying the location of focal 
features, e.g., associated with particular brain functions 
and/or pathologies. Thus topographical and scatter plots 
present different aspects of the fluctuation properties of 
the brain electrical activity of any given subject. Both 
are inefficient for comparison across subjects. What we 
need is a global measure that describes the general, over- 
all feature of all Ui pairs in the form of a single parameter. 

IV. MOMENTS OF THE SCALING EXPONENTS 

We propose to consider the moments of the scaling 
exponents. In general, if we have N numbers, Zj,j = 
1 , • • • , A^, we can calculate the moments 

1 ^ 

where 5 is a positive integer |ic| ]. The information con- 
tained in the first N moments (i.e., q = 1,---,N) is 
enough to reproduce all the Zj by inversion. However, 
we may be interested in only a few of the Gq with lower 
order q, each of which contains some information of all 
the Zj. In our present problem we have N = 128, and we 
shall consider the first ten orders, 1 < g < 10. That is a 
significant reduction of numbers, a process worth inves- 
tigating. 

Before calculating the moments of ai, let us see how 
those values are distributed. Let x be either ai or a2- 
Since no value of has been found to exceed 1.5 in 
the subjects we have examined, we consider the interval 
< a: < 1.5. Divide that interval into M equal cells. 
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which for definiteness we take to be M = 150 here. Let 
the cells be labeled by m = 1, ■ • ■ , M, each having the 
size 5x = 1.5/M. Denote the number of channels whose 
X values are in the mth cell by 7i,„. Define 

Pm^Ura/N. (6) 

It is the fraction of channels whose x values are in the 
range (m ~ l)5x < x < mdx. By definition, we have 
Sm=i Pm = 1. In Fig. 5 we show as an illustration the 
two graphs of P„i for subject A. The two graphs corre- 
spond to ai and a2, and are, in essence, the projections 
of the scatter plot in Fig. 3 onto the ai and a2 axes. 
From Fig. 5 we see that ai is widely distributed, while 
Q!2 is not. Gi gives the average, and G2 is related to the 
width. 

Since the fluctuation of m in should be measured 
relative to its mean, let us consider the normalized mo- 
ments ||ll[ 

m— 1 / \m— 1 / 

(7) 

where i = 1 or 2. Since these moments are averages of 
(m/m)'?, where m is the average-m, they are not very 
sensitive to m itself. They contain the essence of the 
fluctuation properties of ai^2 in all channels. In terms 
of the scaling exponents explicitly, let us use ai{j) to 
denote the value of for channel j so that Eq. (|^) may 
be rewritten as 

In principle, it is possible to examine also the moments 
for q < 0, which would reveal the properties of P„i at 
low values of m. However, the accuracy of our data is 
not too reliable for low-fc analysis, since the 60 Hz noise 
due to ambient electric and magnetic fields has not been 
cleanly filtered out. In this paper, therefore, we restrict 
our study to only the positive q values. For high q, the 

large m/m parts of Pm'^^ dominate Mq^''^\ 

(1 2) 

In Fig. 6 the q-dependences of InMq ' ' are shown for 
the distributions exhibited in Fig. 5 for 2 < g < 10. They 
are approximately linear except for the low values of q. 
The same type of dependencies on q are found for all 
subjects. In Fig. 6 we show two straight lines that can fit 

(i) 

very well the nearly linear behaviors of In Mq vs q for 
g > 5. Thus for large q we have 

M« ocexp q>5. (9) 

The linear extrapolations of the lines to lower values of 
q show the degree of deviation of the the calculated val- 
ues of InMg^'^^ from linearity. Since lnA'4^' and InMf ^ 



behave so similarly in their departures from the linear 
dependencies on g, we plot InM^^^ vs InM^"'^'' in order to 
exhibit their direct relationship without explicit depen- 
dence on q. We find that they are linearly related over a 
wider range of values. This linearity is found to be true 
for all subjects. The plots for three of them are illus- 
trated in Fig. 7, where the straight lines are the linear 
fits. Thus the implication is that there exists a universal 
power-law behavior 

Mf)oc(MW)'' (10) 

valid for all subjects examined. From Eqs.(||) and ( |l0| ) 
we obtain 

f? = M2/Mi, (11) 

but now 7/ is meaningful for all q (except for the lowest 
points) and in that sense independent of q. Thus we 
have discovered a global measure 77 that characterizes all 
ai values of a subject, and varies from subject to subject. 
We postpone the display of all the 77 values for all subjects 
until later. 

To understand the exponential behavior in Eq.(^, we 
note that Gq is dominated by large Zj when q is large, as 
is self-evident in Eq.(^). For asymptotically large q, we 
have Gq oc exp (glnzmax), where ^max = max{zj}. For 
intermediate q, all large values of Zj can make important 
contributions, and the exponential dependence on q can 
still prevail. The denominator in Eq.(0) is Gl = z'^ ~ 
exp ( q Inz ) , where z is the average of Zj , so it is also 
exponential for any q. It is therefore clear that Eq.(H) 
follows, and that fXi depends on all Zj with more weight 
on the large zj values. The power-law behavior of Eq.dlQ) 
implies that the exponent 77 is independent of q and that 
all ai values are relevant contributors to the universal 
behavior. This is an important point worth emphasizing: 
the independence of ry on g implies that the whole spectra 
of ai and 02 are summarized by the one index rj. The 
fact that 77 varies from subject to subject is a consequence 
of the variability of all 128 pairs of (01,012) across the 
subjects, and offers the possibility that 77 can be used as 
a discriminating representation of the brain state. 

V. CORRELATIONS OF THE SCALING 
EXPONENTS 

The analysis in the preceding section treats the mo- 
ments of oi and 02 separately. Only in the last step are 
the global properties embodied in Mg^-* and Mq^^ related 
through the exponent rj in Eq. (|lO|) . In that approach the 
pairing of ai with 02 in each channel is not taken into 
account. However, we know that there are channels, such 
as Ch. 1 in Figs. 1 and 2, where the absence of a dom- 
inant mode of oscillation results in ai « 02. Thus the 
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correlation between the two scaling exponents is an im- 
portant feature that should be explored and quantified. 
To that end we define 

/3 = a2/ai (12) 

for each channel. In most cases we have /? < 1, but (3 > 1 
is possible and, by its rarity, noteworthy. 

Prom a scatter plot, such as Fig. 3, it is possible to 
visualize the /3 distribution, since f3 is just the slope of 
a line from the origin to each point. We show in Fig. 8 
the (3 distributions for the same three subjects as those 
in Fig. 7. Subject B is chosen for display because it has 
the largest 77, while subject C is chosen because it has 
several (3 values that exceed 1. 

To summarize the 128 values of (3j for each subject, we 
apply to them the moment analysis that is developed in 
Sec. IV. Let us therefore define 

The q dependence of In A^^ for the same three subjects are 
shown in Fig. 9. Again, linear fits are very good. Thus 
we have 

Ng (xcicp{ i/q) (14) 

with a distinct u for each subject. Clearly, the ones with 
wide distributions relative to their means have higher 
values of v. 

We now have found two indices, r) and ly, for each sub- 
ject. They describe different aspects of the scaling expo- 
nents. To display those values, it is illuminating to show 
the scatter plot of (ry, i^), as in Fig. 10, which has 18 points 
for the 18 subjects studied. The subjects A, B and C are 
denoted by distinct symbols, same as in Fig. 9. We see 
that the points in Fig. 10 appear to form a band, roughly 
correlated in rj and v. Only subject C has a value of 
that lies above the band, and it is C who has several /3 val- 
ues above 1. Whether that is an anomaly carrying some 
physiological significance is an issue outside the scope of 
this paper, since we assume no knowledge of the physical 
conditions of the subjects. We are similarly ignorant at 
this point about the meaning of the spread of the (ry, p) 
values. However, it is encouraging that the scatter plot 
in (77, i^) is widely distributed for the 18 subjects, since it 
offers the possibility of being a discriminating tool, quite 
different from the alternative scenario, if the parameters 
determined in an analysis had turned out to be nearly 
the same for all subjects. 

VI. CONCLUSION 

Recognizing that the brain is a highly nonlinear sys- 
tem, we have explored a possible way of analyzing the 
EEG time series that avoids the assumptions made in 



linear analyses and in chaos studies. By studying the 
fluctuations from linear trends defined over varying time 
scale, we have found two scaling regions in which the 
RMS fluctuations can be characterized by two dimen- 
sionless scaling exponents, ai and a2, for each channel. 
We then performed moment analyses to reduce the large 
number of pairs of {ai,a2) to simple and manageable 

(i) 

forms. The two types of independent moments, Mq 
and Ng, yield two indices, rj and i/, which provide con- 
cise signatures of the nonlinear behavior of all channels 
of the EEG signals. 

Our emphasis in this paper has been on the nic;thod of 
analysis rather than on the physiological interpretation 
of the results. For the latter task it is necessary to have 
not only more data, but also detailed information on the 
physical conditions of the subjects so as to have a refer- 
ence frame to calibrate the indices obtained. However, 
working with 18 subjects is sufficient to demonstrate the 
effectiveness of the method, to show the universality of 
the scaling behaviors, to reveal the range of variability of 
the indices derived, and to offer the possibility of a new 
way of understanding human brain activity. 
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FIG. 2. F{k) vs k for the three channels in Fig. 1. The 
vertical scales of Ch. 1 and Ch. 2 are shifted upwards by 1.0 
and 0.5 units, respectively. 
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FIG. 3. Scatter plot of a2 vs qi for subject A. The three 
channels exhibited in Figs. 1 and 2 are shown as circles. 
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FIG. 1. A sample of EEG time series in three channels. 
The vertical scales of Ch. 1 and Ch. 2 are shifted upward by 
60 and 30 nV, respectively. 
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FIG. 4. Topographical plots of ai (top) and Q2 (bottom). 
In each figure, ten contour lines are drawn within the data 
range: solid lines above the mean, dotted lines below. 
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FIG. 6. The q dependence of In M^'' for subject A. The 
straight lines are linear fits of the points for q > 5. 
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FIG. 7. A log-log plot of a4^^ vs M^^' for three subjects 
A, B, and C. The solid lines have the slopes given by Eq. (pT|). 
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FIG. 10. Scatter plot of v vs q for all 18 subjects, three of 
which have individual symbols: triangle (A), square (B), and 
star (C). 
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